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Abstract 

We study the velocity of bubble walls in the electroweak phase 
transition. For several extensions of the Standard Model, we estimate 
the friction and calculate the wall velocity, taking into account the 
hydrodynamics. We find that deflagrations are generally more likely 
than detonations. Nevertheless, for models with extra bosons, which 
give a strongly first-order phase transition, the deflagration velocity 
is in general quite high, 0.1 < v w < 0.6. Therefore, such phase transi- 
tions may produce an important signal of gravitational waves. On the 
other hand, models with extra fermions which are strongly coupled to 
the Higgs boson may provide a strongly first-order phase transition 
and small velocities, 10~ 2 < v w < 10 _1 , as required by electroweak 
baryogenesis. 

1 Introduction 

The electroweak phase transition may give rise to a variety of cosmological 
relics such as the baryon asymmetry of the universe, cosmic magnetic fields, 
inhomogeneities or gravitational waves. To be observable, some of these relics 
depend on the strength of the phase transition. This is the case, e.g., of the 
baryon asymmetry of the Universe (BAU) and gravitational waves (GWs). 
Both can be generated in the electroweak phase transition, and both re- 
quire a strongly first-order phase transition. In a first-order phase transition, 
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bubbles of the stable broken-symmetry phase nucleate at a temperature T n 
below the critical temperature T c , and grow inside the supercooled symmet- 
ric phase. The expansion of bubbles provides the required departure from 
thermal equilibrium. 

Gravitational waves are generated by the collisions of bubble walls and 
by the turbulence they produce [U E]- As a consequence, bubble walls with 
higher velocities will generate GWs of larger amplitudes. For electroweak 
baryogenesis [3] , the violation of baryon number must be suppressed in the 
broken-symmetry phase in order to avoid the wash-out of the generated BAU. 
This puts a condition on the vacuum expectation value (VEV) of the Higgs 
field in that phase, namely, (<f>)/T > 1. This quantity is the order pa- 
rameter of the phase transition, and the above condition implies that the 
phase transition must be strongly first-order. The amplitude of the baryon 
asymmetry also depends on the wall velocity. However, unlike the case of 
gravitational waves, the BAU has a maximum for wall velocities in the range 
v w ~ 10~ 2 — 10 _1 , depending on the model [U [5]. For higher velocities the 
sphaleron processes do not have enough time to produce baryons, whereas 
for lower velocities equilibrium is restored and the generated baryon number 
is washed-out. The low velocities required for baryogenesis seem to be too 
small to generate observable gravitational radiation. 

For baryogenesis calculations the wall velocity is often assumed to be of 
the form [S El El El El ED] 

v w « Ap(T)/ V , (1) 

where Ap is the net pressure acting on the wall, i.e., the pressure difference 
between the two phases, and rj is a friction coefficient. This approximation 
corresponds to the case of subsonic walls, which propagate as deflagrations. 
In the context of gravitational waves, on the contrary, the wall is generally 
assumed to propagate as a detonation 0, [TT1 [T21 EE! EH [15J, [El [T7] . Further- 
more, a Jouguet detonation is assumed, leading to a simple expression for 
the velocity [T8] . 

y/I73 + y/a 2 + 2a/3 
Vw ~ l + a ' [) 

which depends only on the parameter a = L/p ra d, where L is the latent heat 
and p ra d is the energy density of radiation. Equation (j2J) does not depend 
on the friction, and the Jouguet velocity gives in fact a lower bound for the 
detonation velocity [121 ED] • In a recent paper [2TJ we investigated the wall 
velocity as a function of L, r], and the supercooling temperature T n . We 
found that Eq. (JTJ is a good approximation for nonrelativistic velocities, 
while Eq. (J2J grossly underestimates the actual detonation velocity. 
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The calculation of the friction coefficient 77 is involved. It depends on 
all the particle species that are present in the plasma and their interactions. 
An accurate computation of rj thus requires considering the details of the 
particle content in each specific model (see e.g. [3 [7]). To compare different 
models it is more workable to consider approximations which depend on a 
few parameters. Such a simple approximation was obtained in Ref. [22] using 
previous results 01 El d El E] , 




for particles with a thermal distribution, whereas infrared bosons contribute 
a term [TU] 

E f^logKW^)- (4) 

bosons 

In these equations, gi is the number of degrees of freedom (d.o.f.) of species 
i with Higgs-dependent mass = hi<f), Tj are interaction rates which are 
typically < 10 _1 T, Xi — 2 for fermions and Xi — m i (0) /T f° r bosons, a is 
the surface tension of the bubble wall, m 2 D ~ h]T 2 is the Debye mass squared, 
and L w is the width of the bubble wall, L w « <p 2 /&- The derivation of Eqs. 
(J3j) and (TJ1) involves expanding the distribution functions to lowest order in 
rrii/T, and thus they break down for rrii/T 1. In particular, for a squared 
mass of the form m 2 — fi 2 + h 2 <f) 2 , with large /ij, the particle density in the 
symmetric phase will be suppressed by a Boltzmann factor exp(— /ij/T), and 
this species will not contribute to the friction. 

In this paper we study the velocity of bubble walls in the electroweak 
phase transition for several extensions of the Standard Model (SM). For that 
aim, we modify Eqs. and (jlj) to take into account more general masses 
mj (</>). We also take into account the effects of hydrodynamics in order to 
consider both deflagration and detonation solutions. We consider several ex- 
tensions of the SM, including the Minimal Supersymmetric Standard Model 
and extensions with singlet scalars and with heavy fermions. We also inves- 
tigate the effect of cubic terms in the tree level potential. 

The plan is the following. In section [2] we review the dynamics of a first- 
order electroweak phase transition. In section [3] we find an approximation 
for the friction coefficient which is valid for large as well as for small values of 
rrii((j))/T. In sectionHJwe write down the equations for the wall velocity, which 
we solve numerically. The result depends on several parameters, namely, 
the critical temperature T c , the nucleation temperature T n , the latent heat 
L, and the friction coefficient rj. We compute these parameters for several 
models in section [3, and we calculate the wall velocity. Finally, in section [6] 
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we discuss the implications of our results for baryogenesis and gravitational 
wave production in the electroweak phase transition. Our conclusions are 
summarized in section [7J 



2 Dynamics of the electroweak phase transi- 
tion 

In the SM, the electroweak phase transition is only a smooth crossover. How- 
ever, many extensions of the model give a first-order phase transition. For 
simplicity we shall consider models with a single Higgs field, or models for 
which considering a single Higgs provides a good approximation. Thus, our 
theory will consist of a tree-level potential 



V 



-m 2 (b 2 + - ( 
4 



(5) 



for a scalar field (the background Higgs field, defined by (H°) = 0/\/2). 
The vacuum expectation value of the Higgs is given by v — a/2/ Am = 
246GeV, and A fixes the Higgs mass, m 2 H = 2Xv 2 . Imposing the renormal- 
ization conditions that the minimum of the potential and the mass of do 
not change with respect to their tree-level values [23], the one-loop zero- 
temperature potential is given by V (0) = Vq (0) + V\ (0), with 



Vi (0) = £ 



±9i 
64tt 2 



mf{(f)) lo 



m, 2 (0) 



2m-(0)m-(u) 



+ c, (6) 



where gi is the number of d.o.f. of each particle species, m« (0) is the 0- 
dependent mass, and the upper and lower signs correspond to bosons and 
fermions, respectively. We have added a constant c such that V (v) = 0, so 
that the energy density vanishes in the true vacuum at zero temperature. In 
the symmetric phase we will have a false vacuum energy density, given by 
Pa = V (0), which contributes to the Hubble rate during the phase transition. 
For particle masses of the form mf = h 2 (j) 2 + /j, 2 , we have 



Pa 



T9i 
32tt 5 



h 



2/z? (* 



2i tH 

v 



V 



T (7 > 



The free energy is given by the finite-temperature effective potential. To 
one-loop order, including the resummed daisy diagrams, we have 



^ = y o (0) + F 1 (0)+ ^(0,T), 
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where the finite-temperature corrections are given by 



' 2?r 2 
ft?" 



(ix a; 2 log 



1 =F exp -\ x 2 + m 2 (0) /T : 



bosons 



(9) 



where the upper sign stands for bosons, the lower sign stands for fermions, 
and M 2 ((f)) = m 2 (0) + ITj (T), where II j (T) are the thermal masses. The last 
term receives contributions from all the bosonic species except the transverse 
polarizations of the gauge bosons. 

At high temperature the symmetry is restored, and in a certain range 
of temperatures, the symmetric minimum = coexists with a symmetry- 
breaking minimum <j) m (T). The free energy density of the unbroken-symmetry 
phase is given by J- U (T) = JF(0,T), whereas that of the broken-symmetry 
phase is given by TiiT) = jF(0 m (T),T). The critical temperature is that 
for which TuiT^) = JF 6 (T C ). The energy density in each phase is given by 
p (T) = T(T) — TJ-'(T) (a prime indicates derivative of a function with 
respect to its variable), and the latent heat L = p u (T c ) — p b (T c ) is given by 



L = T c {F b {T c )-F' u {Tc)). 



(10) 



We define the thermal energy density p u by subtracting the vacuum energy 
density, p u = p u — Pa- In general, we have p u ~ 7r 2 g*T 4 /30, where is 
the number of relativistic d.o.f. The hydrodynamics of the bubble wall will 
depend on the parameters a c = L/p u (T c ) and a n = L/p u (T n ) , where T n is 
the nucleation temperature. 

The nucleation of bubbles [23, 26J is governed by the three-dimensional 
instanton action 



S, =4tt 



r dr 



\\f r ) 



where 

V T {4>)=H^T)-T{Q : T). (12) 

The bounce solution of this action, which is obtained by extremizing S3, gives 
the radial configuration of the nucleated bubble, assumed to be spherically 
symmetric. The action of the bounce coincides with the free energy of a 
critical bubble (i.e., a bubble in unstable equilibrium between expansion and 
contraction). This solution obeys the equation 



cP 



+ 



dr 2 r dr 



(13) 
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with boundary conditions 

^r(O) = 0, lim (r) =0. (14) 

We will solve Eq. (jT3l) iteratively by the overshoot-undershoot method]!]- The 
thermal tunneling probability for bubble nucleation per unit volume and time 
is [26] 

r n (T)^A(T)e^ T » T , (15) 

with A (T) = [S3 (T) / (2nT)] 3 ^ 2 . The nucleation time t n is defined as that at 
which the probability of finding a bubble in a causal volume is 1, 

dtT n (T)V c = l, (16) 



where t c is the time at which the Universe reached the critical temperature T c 
and, in the radiation-dominated era, the causal volume is given by V c ~ (2t) 3 . 
The time-temperature relation is given by 

dT/dt = -HT, (17) 



where H is the expansion rate, H = ^d>iiGp u (T) /3. Here, G is Newton's 
constant. Using Eq. (|T7j) we can solve Eq. ({TBI) for the temperature T n 
at which the first bubbles are nucleated. If p u ~ p u »s 7r 2 tyr =|t T 4 /30, then 
the time-temperature relation is given by the usual expression t = £Mp/T 2 , 
where Mp is the Planck mass and f = ^45/(16tt 3 #*). 

It is useful to consider the profile (f> (r) of the critical bubble at T w T c . 
At the critical temperature, the radius of the nucleated bubble diverges (and 
the nucleation rate vanishes). Hence, for T w T c the second term in Eq. f|T5|) 
can be neglected, since the wall is much thinner than the radius. Thus, one 

obtains 

d(j>/dr = -^2V T {<t>{r)). (18) 

Within this approximation, the wall is planar and its surface tension a = 
J-™ {d(p/dr) 2 dr is given by 



a 



(19) 



where <p c = m (T c ). We can also invert relation ( |T8l) to obtain the wall width 
L TO . If we define, e.g., L w = r (0 = O.10 c ) — r (0 = O.90 c ), we have 

/■0.9(/)c 

L w = / d((>/y/2VT ((/>). (20) 



L See Ref. 27 for details. 
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Roughly, a and L w are related by a ~ <pl/L w , and the wall width is given by 

L tt ~y/(o)- 1/2 . 

3 Microphysics 

According to kinetic theory, for a planar wall in stationary motion along the 
z direction, the friction force per unit area is given by [6] 

fa* on = £ ft y ^ y _3 ___*/„ (2D 

where -E 1 = \/ p 2 + m 2 and 5/j is the departure from the equilibrium distri- 
bution f (Ei/T) for each particle species, with 

/°(*) = -TTT- (22) 
e x ± 1 

The kinetic description is valid for particles with p 3> L~ l , for which the 
background field varies slowly and the semi classical (WKB) approximation 
is valid. Since in general L w l T, this condition is satisfied for all but 
the most infrared particles 0|. As usual, we will assume that the friction is 
proportional to the wall velocity, 

friction = r)V w , (23) 



where the friction coefficient r\ is obtained by considering Eq. ( |2Tj) to linear 
order in v w . The deviations 5 fa can in principle be calculated by considering 
the Boltzmann equation for the distribution functions [H El El [HI IH] ■ How- 
ever, infrared excitations of bosonic fields should be treated classically and 
undergo overdamped evolution [TQl EEJ [29]. We shall refer to these fields 
as "infrared bosons", whereas we shall call "thermal particles" those which 
obey the Boltzmann equation. 

3.1 Thermal particles 

We begin by considering the case of thermal particles. It is usual to em- 
ploy the ansatz f = fo (E/T — fi/T + E5T/T 2 + p z v/T) for the distribution 
functions [HI U\- In that case, a system of equations for fi, 5T and v for each 
particle species can be derived from the Boltzmann equation 

[d t + (d p E) d z - (d z E) d Pz ] f = -C [f] , (24) 
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where C [/] is the collision term. We need to simplify further the problem 
in order to obtain a simple analytical expression which can be applied to 
different models. Therefore, we shall use the ansatz f = fo (E/T — 5), which 
is equivalent to considering only the term fi/T. Hence, the deviation from 
fo(E/T) is 5f = —^(E/T) 5, and the equation for 5 is obtained by lin- 
earizing the Boltzmann equation. Assuming stationary motion and making 
a momentum integration, one obtains, for each particle species (see Ref. [22] 
for details), 

d5 civ w dm 2 
C2 ^-™ = 2T^7' (25) 
where T is an interaction rate arising from the collision integrals, and the 
coefficients c\ and c<i are defined by 



It is out of the scope of this work to calculate the collision integrals, which 
depend on the particle content of each model. Numerically, the rates V are 
~ l(r 2 T [6], and we shall set Y/T = 5 x KT 2 . In the thick wall limit, the 
first term in Eq. ( 1251) can be neglected, and we obtain 

civ w dm 2 ^ 



2T 2 r dz 



We now insert 5f = -f' Q (E/T)8 in Eq. ([21]), with 5 given by Eq. (J2j 
with m? = fi 2 + h 2 (p 2 for each particle species i. Performing the momentum 
integration we obtain 

l4 r+oo 
xvH I 2 i2 1 12 



V = 2^-jr I tiWdz. (28) 



3C 



For T w T c we can use the thin wall approximation^! <$' = \/2Vt ((f)) in Eq. 
(I28p . Thus, the friction caused by particles with thermal distributions is 
given by 



= l^ 9J ¥ \ ' ZMtfy/WM- (29) 



Particles with larger couplings h give the main contributions to the friction 
force, since they have stronger interactions with the bubble wall. 



2 The wall is thick in comparison with T , but it is thin in comparison with the bubble 
radius. 
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3.2 Infrared bosons 



Infrared boson excitations must be treated classically [28] and undergo over- 
damped evolution |29j. Relating the population function to the squared 
amplitude of the field [TO], one obtains the equation 

^4 = -EH f , (30) 
8p at 

where m? D is the squared Debye mass, m 2 D = (ll/6)g 2 T 2 for the W and 
Z bosons of the SM, and m 2 D = h 2 T 2 /3 for a scalar singlet. Writing / = 
fo (E/T) + Sf we have, to first-order in v w , 

r r Km 2 D dm 2 

Sf= -m^ fo ~d^ Vv " (3) 

Inserting the departure from equilibrium f[3"TT) in Eq. f[2~Tl) and doing the 
momentum integration, we obtain 

1= £ 9^h* [ +0 ° b^dz. (32) 

bosons 

with 

6 = _ T 2 /" f'JE/T). (33) 

7 (2tt) 3 pE 4 ov 7 ; v ; 

The integral in Eq. (1321) will diverge if m(4>) vanishes. Indeed, for small 
m/T, the momentum integral (1331) is infrared dominated. Therefore, we can 
make the approximation f (x) ~ —1/x 2 , and we obtain 

b = T 4 / (87r 2 m 4 ) . (34) 

Hence, for \i = the integral in Eq. (I32p has a logarithmic divergence at 
= 0. However, the kinetic description (12 ip breaks down for very infrared 
particles, which in the case m — > dominate. The actual contribution of 
degrees of freedom with p < L~ l is subdominant, because their wavelength 
cannot resolve the thickness of the wall (TO]- Thus, for small \x we will cut 
off the integral at m (0) = L~ l , and we have 



E ^-^k h ' dM?VWrb (^) , (35) 

" J (bn 



X — T ; 1 ' ' ! ' 

Vir = 

bosons 



with 0o — V 'L~? — jj? /h for /i < L" 1 , and O = for /i > L^ 1 . For a particle 
with small enough /i, we will have 0o > C - In such a case the particle will 
not contribute to the friction and we shall set r\ — 0. 
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Since Eqs. ( 1291) and (1551) depend on the profile of the wall, the contri- 
bution of a particle species to the friction depends on the whole particle 
content of the model. As an example, we consider the case of the SM with 
an additional complex scalar field. Figure Q] shows the contribution of the 
extra singlet to the friction as a function of the coupling h. The dashed line 
corresponds to the friction coefficient r] th given by Eq. (12U1) . the dotted line 
to ?7i r given by Eq. (135|) . and the solid line shows their sum. In Fig. [5] we 
fix h — 1 and vary [i. The coefficients 77th and 7]- ir dominate in different pa- 
rameter regions, and we shall use rj = rfth + rfa for applications in this paper. 
As expected, rj ir dominates for relatively small values of /i and h. We see 
that the friction increases with h, as particles interact more strongly with 
the wall, and decreases with \i due to Boltzmann suppression of the particle 
density. 
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Figure 1: The contribution of a complex singlet to the friction as a function 
of the coupling h for fi = and win = 125GeV. 



3.3 Limiting cases 

In section [5] we shall compute the friction coefficient for each species numer- 
ically. It is useful, though, to consider here some limiting cases. 

The coefficient c\ depends on m ((f)). It is usual to consider only the lowest 
order in m/T. In that limit we have 

Cl/ = log2/27r 2 , ci6 = -log (m/T) /2vr 2 , (36) 
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Figure 2: The contribution of a complex singlet to the friction as a function 
of n for h = 1 and m# = 125GeV. 

for fermions and bosons respectively. Since cu depends only logarithmically 
on m, we may replace m (0) ~ m (<f) c ) in Eq. f[3"oT) . The integral in Eq. f[2"9"]) 
goes like 2 cr, and in the small m/T limit we recover Eq. (jHJ), 

^~^&feV for m/T «1, (37) 



r v 2tt 2 ^ 

with x = 2 for fermions and x = \/^ 2 0c + /i 2 /^ ~ ft- for bosons. However, 
the approximation (IHUj) will break down whenever /x or /i0 become large in 
comparison with the temperature. In particular, particles with large h in 
a strong phase transition will have m ~ /i</> ^> T in the broken-symmetry 
phase. 

In the limit m/T ^> 1 we have 



ci 



(m/T) 1/2 exp {-m/T) / (2tt) 3/2 . (3* 



It is apparent that the contribution of a heavy particle to the friction is 
suppressed by a Boltzmann factor. In the case of small \x but h<p c ^> T, we 
can estimate the integral in Eq. (T2"9"l) using a quartic approximation for the 
effective potential at T = T c , V T (0) « A0 2 (0 C - 0) 2 with A = V£ (0) /(20 2 ). 
Since for small the integrand is suppressed by a factor 3 , we can use the 
approximation ( 138]) in Eq. ( 1291) . To leading order in h(p c /T c , we obtain 

Vth « p for /i0 c /T c > 1. (39) 
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We see that in the limit of a very strong phase transition the friction coeffi- 
cient does not increase like h A as in Eq. (|37[) . This is because the integrand 
in Eq. ( l29l) is suppressed for hep > T. Still, r\ grows with the strength of 
the phase transition due to the factor sJV^ (0). For a heavy particle with 
H ^> T, the function c\ does not depend on 0, and the integral in Eq. ( 1291) 
goes with c\cj) 2 c <j. Therefore, we have 



9h * ^ e -^ T for [i I T » 1. (40) 



(2tt)' r T 



This exponential suppression of the contribution to the friction is important, 
since the heavy particle will also be decoupled from the effective potential. 
For instance, as the value of [i for a boson is increased, the phase transi- 
tion becomes weaker. Consequently, there will be less supercooling and a 
smaller pressure difference at T = T n . Nevertheless, the wall velocity will 
not necessarily decrease, since the friction will also be smaller. 

For the coefficient b, we have the approximation (1341) for m/T <C 1. For 
small we also have Vr (0) ~ (0) 2 , and we obtain 



V ~ 64tt 



' m(0 c ) (m 2 (0 c ) -m 2 (0 o ))/i 2 
log 



for m/T <C 1. 



m (0 O ) 2m 2 (0 C ) m 2 (0 O ) 

(41) 

The factor of h A in (1351) has disappeared, but a factor /i 2 still remains in 
m 2 D . This is the generalization of Eq. (j4j) for small but nonvanishing \i. 
Notice that for \i < L~ l we have m(0 o ) = L" 1 , while for > L" 1 we have 
m(0 o ) = A*- Thus, noting that \/Vfj0) ~ L" 1 , we recover Eq. (T4]) in the 
limit \i — > 0. On the other hand, for /i 7^ the friction is smaller. 
In the limit of very large m/T, we have 

6 = e - m / T T 3 /(27T 2 m 3 ). (42) 

Let us consider, e.g., the case /i = and h(p c ^> T. Since 6 is exponentially 
suppressed for m = h<f) ^> T, we cut off the integral in Eq. (I35I) at = 
T/h <p c and use the small m approximation for 6. Thus, we obtain 

gml^V£ (0)T C 

r? 7:1 log(TL^) for ft0 c /T c > 1. (43) 

047T 

Due to the infrared behavior, the result is similar to the case /1 = 0, hcf) c /T c <C 
1. The only difference is that the log is evaluated at T instead of m(0 c ). For 
a particle with large [i the contribution is 

_ gh'ml^/T for/VT>>1; (44) 



167T /i 



3 
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which is exponentially suppressed, as expected. 

In Figure [3] we have plotted the total friction coefficient (including the 
contributions of the top quark and gauge bosons) for the case of the SM with 
an additional singlet, as a function of the coupling h of the extra singlet. 
The dashed line corresponds to using the approximations given by Eqs. ([3]) 
and (j3J). Due to numerical factors in these rough approximations, this curve 
does not match the solid line for small h. Apart from this fact, the behavior 
is similar. However, for larger vales of h the dashed line grows faster and 
crosses the solid line. This is because the approximation increases with h A , 
while the correct result does not, according to Eq. (1391) . 



10° 
10" 1 



10~ 3 
10" 4 

0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 

h 

Figure 3: The friction coefficient for the SM with a complex singlet as a 
function of the coupling h for m H = 125GeV and /x = is plotted in solid 
line. The dashed line corresponds to using the small m/T approximations. 

We wish to emphasize the fact that the friction depends on the wall profile 
besides depending on the parameters g and h. For instance, in the case of 
small m/T, r\ is proportional to (j) 2 c a ~ (p^L' 1 . Sometimes the wall width 
is estimated as L w « 1/T. However, L w may change significantly with the 
strength of the phase transition. We have computed L w and a numerically 
using Eqs. ( |T9l) and ( 1201) . For the parameters used in Fig. [3] they vary 
from L w ~ 100/T and o ~ 10~ 4 T 3 for small h (weak phase transitions) to 
L w ~ 1/T and a ~ T 3 for large h (strong phase transitions). This dependence 
will become particularly important if the strength of the phase transition is 
varied without changing the particle content of the model (e.g., by changing 
the tree-level potential). 
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4 The bubble wall velocity 



If we ignore hydrodynamics, the steady state velocity of the bubble wall 
is obtained by equating the friction force (1231) to the pressure difference 
between phases. We shall use subindexes b and u for quantities in the 
broken and unbroken symmetry phase, respectively. Thus, ignoring hy- 
drodynamics, we have r)V w = Pb{T) — p u (T). Here, the pressure difference 
Ap(T) = pb (T) — p u (T) is given by the effective potential ( |T2l) . Ap(T) = 
—V (4> m {T) , T). However, the motion of the wall in the plasma causes vari- 
ations of the temperature and velocity of the fluid, and is in turn affected by 
these perturbations. 

For hydrodynamic considerations (see, e.g., [3UJ [ITI] EH] ) we can assume a 
thin wall, and the temperature and fluid velocity turn out to be discontinuous 
at the interface. Thus, we have different temperatures T u and T& on each 
side of the wall. We assume a planar wall in stationary motion in the z 
direction. In the rest frame of the wall, the continuity conditions for energy 
and momentum fluxes give the relations [3T] 

Wul 2 u v u = WbJ%v b , (45) 

Wullvl+Pu = Wb^vl+Pb, 

where v is the fluid velocity 7 = 1/ yl — v 2 , and w = p + p is the enthalpy 
density. In Eq. (|4"5"]) we have, e.g., p u = p u (T u ) (notice that the thermo- 
dynamical quantities have in general different values in each phase, even for 
T u = Tb). To obtain a macroscopic equation involving the friction, one can 
introduce a damping term of the form u^d^cf) in the equation of motion for 
the Higgs field. This equation must be integrated taking into account the 
temperature variation. Assuming a thin wall, one obtains [21J 

1 71 

Pu - Pb - 2 ( s « + s b) ( T u ~T b ) + - (v u -f u + v b %) = 0, (46) 

where s is the entropy density and 77 is the friction coefficient obtained from 
the microphysics calculation, r] = [Ap(T) / v w ] micvo . The thermodynamical 
variables are related by an equation of state (EOS), so Eqs. (jl5"|) and (pffi"]) 
have only four unknowns, namely, the velocities v u ^ and the temperatures 
T U £. Besides, the temperature T u outside the bubble can be determined by 
computing the nucleation temperature. 

The velocities v u and Vb of the fluid are measured in the reference frame of 
the wall. Equivalently, |i> u | and \vb\ give the velocity of the wall with respect 
to the fluid on each side of the wall. We want to calculate the wall velocity 
v w in the reference frame in which the fluid is at rest very far in front of the 
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wall and very far behind the wall (at the center of the bubble). We shall 
refer to this reference frame as the laboratory frame. 

The stationary motion of the wall admits two kinds of solutions, called 
detonations and deflagrations. For detonations we have \vb\ < \v u \ and 
| v u | > c s , where c s = l/v3 is the speed of sound in the relativistic fluid. 
Therefore, in this case the wall moves supersonically with respect to the 
fluid in front of it. Hence, no information on the motion of the wall is trans- 
mitted into this fluid, which can be assumed to be at rest in the laboratory 
frame. According to this boundary condition, the detonation is supersonic, 
v w = —Vu > c s , and the temperature T u is the nucleation temperature T n . 
In the reference frame of the wall, the incoming flow is supersonic. The out- 
going flow inside the bubble has a lower velocity and could in principle be 
subsonic (strong detonation). This possibility, however, is forbidden by the 
boundary condition at the center of the bubble. Therefore a detonation can 
only be weak > c s ) or Jouguet (\vb\ = c s ). In general, though, we will 
only have weak detonations [2T] . 

For deflagrations we have \v\\ > \v u \ and \v u \ < c s , so the wall is subsonic 
with respect to the fluid in front of it. It turns out that a single front is 
not enough to satisfy the boundary conditions in this case. Thus, the phase- 
transition front must be preceded by a shock front. In the laboratory frame, 
the fluid between the two fronts has a finite velocity, and outside this region 
the fluid is at rest. Therefore, we have v w = —Vb, and the fluid in front of the 
wall moves in the direction of the latter. In the frame of the wall, there is a 
subsonic flow coming from the symmetric phase, and an outgoing flow with a 
larger velocity v w . In principle, the outgoing flow can be supersonic (strong 
deflagration), although in general we will have v w < c s (weak deflagration). 
The case v w = c s is called a Jouguet deflagration. In the deflagration case, 
the temperature T u is higher than T n , since the fluid in the region between 
the bubble wall and the shock front is compressed and reheated. The relation 
between T u and T n can be obtained by considering the fluid conditions f j45l) 
for the shock front [31J. 

In order to solve Eqs. f|4"5l) and fj4~6l) . it is convenient to use the bag 
equation of state 

Pu (T) = a u T 4 + e, p u (T) = a u T 4 /3 - e, (47) 
p b (T) = a b T 4 , p b (T) = a fe T 4 /3, 

which is the simplest EOS that keeps the essential features of a phase transi- 
tion. The coefficients a u and a& must be different, since in this model we have 
db/ a u = 1 — 3a c , with a c = ej (a u T c 4 ). The constant e is the false vacuum 
energy density of the high-temperature phase, and also determines the latent 
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heat L = p u (T c ) - p b (T c ) b}^ 

L = 4e. (48) 

In the general case, the latent heat and the vacuum energy density will not 
have this simple relation. In order to apply the results of this calculation to 
a general model, it is convenient to rewrite the parameter a c as 

where p u (T) = a u T A is the thermal energy density in the high-temperature 
phase. 

Using the EOS (@7D in Eqs. flS} and flU, we eliminate T b , 



It 



a u ( , 1 + v u v b 



a b \ u l/3- v u v b 



-,1/4 

(50) 



where a u = ej (a u T*), and we still have two equations for v u , v b and a u . One 
of them comes from hydrodynamics alone 



_L _|_ Vk _|_ ,/ / J_ _i_ Vh\ 2 _i_ r.2 _i_ 2 _ 1 
5ui. ^2^1/1 6«i. ' 2 / ' u u ' ft"" 3 



6s)ft 2 \/ \ 6vf, 2 . 

.. V : , , (51) 

where the + and - signs in front of the square root correspond to detonations 
and deflagrations, respectively. The second equation involves microphysics 

4v u v b a u 2/ s b \ ( T b \ 2a u r] 

-- H 1 - — H — {\v u \ 7 M + \v b \ 7 6 ) = 0, (52) 



1 - 3w u w b 3 V s u / V T t 

with Sfe/s u = (a b /a u ) (T b /T u ) 3 . In the case of detonations, we can solve these 
equations for | function of a U7 since a u depends on T u = T n . In 

the case of deflagrations, we have an additional equation relating T u and T n , 

V3(«n-««) V u -V b 

(53) 



A/(3a n + a u ) (3a u + a n ) 1 - ^6 



where a n = ej (a u T*), and we can solve for \v b \ = v w . In any case, T u can be 
eliminated, and the result depends on 

L . . 

an = -rz , r . 54 
4p u (T n ) 



3 A simpler model with a single coefficient a u = a b = a is sometimes used. In that case 
we would have L = e. However, such a model is not suitable for a phase transition, since 
the critical temperature [for which p + (T c ) — p_(T c )] does not exist [21j . 



16 



These equations can be solved numerically [21] . and the wall velocity fi- 
nally depends on the parameters a c and a n , and on the ratio r]/L. Depending 
on the parameters, there may be only deflagrations, only detonations, both, 
or none. The deflagration solution always exists if the friction is large enough 
or the supercooling is small enough. The detonation solution will exist if the 
friction is small enough and the supercooling is large enough. In general, 
there is no detonation solution when the deflagration velocity is smaller than 
w 0.5. As the velocity approaches the speed of sound, the detonation solu- 
tion may appear. The deflagration solution disappears soon after becoming 
supersonic. Analytical approximations were also found in Ref. [2T] for both 
solutions. In this paper we shall use the numerical results. 

Notice that the definitions of a c and a n [Eqs. (jUJ) and (j53l) ] involve the 
thermal energy densities at the temperatures T c and T n , but the latent heat is 
in both cases the energy density discontinuity at T = T c . These parameters 
arise as a consequence of the use of the simple bag EOS, and their definitions 
must be respected in applications. As discussed in Ref. [21], if we used for 
a n the energy density that is released at T = T n , which is larger than L, 
we would be overestimating the velocity, since this would be equivalent to 
considering a stronger supercooling [i.e., a smaller value of p u (T n )]. 

5 The electroweak wall velocity 

In this section we calculate the value of the wall velocity in the electroweak 
phase transition. For that aim, we compute the parameters a c , a n , and rj/L 
for several models. The temperature T n is defined through Eq. (TTB1) . and is 
the temperature at which bubbles begin to nucleate. By the time bubbles 
occupy all space, the temperature will be in general different. If the wall 
velocity is relatively large, bubbles will not have time to interact with each 
other until they percolate. In this case the temperature will just decrease 
due to the expansion of the Universe, and the velocity will increase during 
the phase transition. In general, though, the transition will be short enough, 
so that the velocity will not change significantly from v w (T n ) [22]. On the 
other hand, in the case of slow deflagrations, the supersonic shock fronts 
preceding the walls will cause a reheating during bubble expansion. Hence, 
the temperature will grow and the wall velocity will decrease from v w (T n ). 
Depending on the amount of latent heat, the reheating temperature may be 
very close to T c [221 [23 EH]- In this case, the wall velocity will decrease 
significantly. We shall not take into account this possibility in the present 
paper. We remark, though, that it may have important consequences [33|l3l]. 
We shall consider several extensions of the SM. The relevant SM contri- 
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butions to the one-loop effective potential come from the Z and W bosons, 
the top quark, and the Higgs and Goldstone bosons. It is usual to ignore 
the Higgs sector in the one-loop radiative corrections. This should be a good 
approximation in extensions of the SM which include particles with strong 
couplings to 0. The ^-dependent masses of the weak gauge bosons and top 
quark are of the form hi<p, with hi = rrii/v, where are the physical masses 
at zero temperature. We shall ignore, as usual, the longitudinal components 
of the weak gauge bosons, which are screened by plasma effects. Thus, the 
W and Z contribute corrections of the form ( 161191) to the free energy (jSJ), 
with 4 and 2 bosonic d.o.f., respectively. The top contributes with g t = 12 
fermionic d.o.f. The rest of the SM particles have hi <C 1 and only con- 
tribute a ^-independent term — 7r 2 gi ight T 4 /90, with gn ght « 90. To strengthen 
the electroweak phase transition, extra particles are usually added to the SM, 
with strong couplings to 0. 

5.1 SM extensions with scalars 

It is well known that the easiest way of strengthening the phase transition is 
by extending the scalar sector of the SM. The simplest extension consists of 
adding gauge singlet scalars [231 EH] ■ 

5.1.1 A complex singlet 

We shall first consider a model in which a single complex scalar S is added 
to the SM [111 Ell EE]- The coupling to the Higgs is of the form 2h 2 S^SH^H, 
and the field S may have a SU(2) x U(l) invariant mass term fi 2 S^S and 
a quartic term As (S^S) . For simplicity, we will ignore the possibility that 
cubic terms exist in the tree-level potential, and we will assume fi 2 > 0. 
We notice however that a negative value of fi 2 , as well as a cubic coupling, 
may enhance the strength of the phase transition through tree-level effects 
[3H] . We shall study the effect of a cubic term at the end of this section. 
Hence, the field S gives contributions to the free energy of the form ( I6l9i) . 
with g = 2 d.o.f. and a mass m 2 (0) = h 2 (p 2 + fi 2 . The thermal mass is 
given by n = {h 2 + X s ) T 2 /3 [37J. We shall take As = for our numerical 
calculations. We have checked that considering As 7^ does not introduce 
qualitative differences in the results. 

FigureHlshows the value of the wall velocity at the nucleation temperature 
T n , as a function of the coupling h for different values of the parameters. 
We have plotted both the deflagration and the detonation solutions, when 
they exist, for the case \i = (which maximizes the strength of the phase 
transition) and Higgs masses m# = 100, 125, 150 and 175GeV, and for the 
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case /i = lOOGeV and = 125GeV. For the same h, larger values of 
ran give weaker phase transitions and, consequently, lower values of v w . The 
crosses indicate the points where (f) (T n ) /T n = 1. As h is increased, the phase 
transition becomes stronger. As a consequence, the temperature T n decreases 
and the supercooling stage lasts longer. In this model, if h is too large the 
Higgs VEV remains stuck in the symmetric phase, and the Universe enters a 
period of inflation. We are not interested in such extreme cases, and we just 
plot the curves up to a value of h for which the time required to get out of 
the supercooling stage becomes too long for the numerical computation. 




h 

Figure 4: The wall velocity as a function of the coupling h of the singlet, for 
several Higgs masses. All the curves correspond to the case /i = 0, except 
for the doted line, which corresponds to /i = lOOGeV. The crosses indicate 
the case (p n /T n = 1. 

Although the amount of supercooling grows with h, so does the friction 
and the behavior of the wall velocity is rather unpredictable. As can be seen 
in Fig. HI the velocity in general increases with h. However, for larger values 
of m H (e.g., the case m H = 175GeV), the velocity eventually decreases for 
large h. This is because larger values of h are required to achieve a phase 
transition of the same strength, and the friction becomes important. Only 
the cases m# = 100 and 125GeV have detonations (the short curves which lie 
above the speed of sound). Notice that there is a small range of h for which 
there are neither deflagrations nor detonations, i.e., no stationary solution for 
the wall velocity. In such a case, the steady state will not be reached, and the 
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wall will accelerate until the end of the phase transition. The case in which 
the wall velocity gets close to the speed of light may have important implica- 
tions for gravitational wave generation. Such ultrarelativistic velocities have 
been considered recently [39] for the SM extension with singlet scalar fields, 
finding that "runaway" solutions exist for very strong phase transitions. 

For /i 7^ the phase transition is weaker and the amount of supercooling 
is lower, but the friction is lower too. As an example, consider the case 
m H = Y2hGeV and fi = lOOGeV (dotted line in Fig. HJ). We see that for 
small h the velocity is lower than in the case of \x = 0, while for large h the 
velocity is higher. In Fig. [5] we have plotted the wall velocity as a function 
of \i for run = !25GeV and h = 1.4, 1.6 and 1.8. We see that for small \i the 
wall velocity increases with /i, indicating that the friction decreases faster 
than the strength of the phase transition. For large \x the velocity eventually 
decreases, as the extra boson decouples from the thermal plasma. For the 
cases h = 1.4 and h — 1.6 there are only deflagrations. The crosses in the 
curves indicate a phase transition with (T n ) /T n = 1. To the right of this 
point the phase transition is weaker. The value h = 1.8 is quite large, and 
there are detonations for a long range of values of /i. In this case, deflagrations 
appear only for fi > 170GeV. We see that there is a range of parameters for 
which there exist both kinds of solution for the wall propagation, and also a 
range in which there is no solution. 
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Figure 5: The wall velocity as a function of fj,, for m# = 125GeV and three 
values of h. 
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5.1.2 A hidden sector 

Recently, an extension of the SM with several real singlets Si has been con- 
sidered [40J. These bosons constitute a hidden sector which couples only to 
the SM Higgs doublet through a term h 2 H^H ^ Sf (for simplicity, universal 
couplings hi = h are assumed). Following Ref. [ID], we assume there are no 
linear, cubic, or quartic terms in the hidden-sector scalar fields. Therefore, 
this case is similar to the previous one, with m 2 (0) = h 2 (p 2 + ji 2 and g d.o.f., 
where g is the number of singlets. We shall take g = 12, as in Ref. [40J. 
If the fields Si do not have mass terms, so that they only get a mass from 
electroweak breaking, the phase transition can be made exceedingly strong. 
Interestingly, this model allows to consider the classically conformal case, 
corresponding to m 2 = in Eq. (jSJ). The loop corrections break conformal 
invariance, and a mass scale appears via dimensional transmutation. Im- 
posing appropriate renormalization conditions, for a given Higgs mass the 
classically conformal case occurs for a fixed value of h. Since we are inter- 
ested in the variation of the strength of the phase transition as a function of 
the parameters, we shall not give special attention to this particular case. 

Figure [6] shows the velocity v w (T n ) as a function of the coupling h for 
jj, — and different values of m# in the range 100 — 2QQGeV . Since this model 
has more d.o.f. than the previous one, the phase transition becomes strongly 
first-order for lower values of h. As a consequence, the wall velocity grows 
more quickly with h. However, we see that in all the cases the nucleation 
stage becomes too long before detonations can exist. For /i ^ the phase 
transition weakens and higher values of h can be considered. We have found 
that detonations appear in the range [i ~ 50 — X'hQGeV for low values of mn 
and extreme values of h. For most values of h there are only deflagrations, 
and the behavior with \x is similar to that of the lower curves of Fig. El 

5.2 The MSSM 

The Minimal Supersymmetric Standard Model (MSSM) has been extensively 
investigated in connection to electroweak baryogenesis [H], since it may pro- 
vide a strong phase transition and sources of CP violation. The MSSM 
contains two complex Higgs doublets Hi and H 2 . We define the vacuum 
expectation values v\ = (H®) and v 2 = (H®)- It is customary to simplify 
the problem by considering the limit in which the CP-odd Higgs mass is 
large (m^ ^> rriz). In this limit the low energy theory contains a single 
Higgs doublet $, and the masses and couplings depend on tan/3 = v 2 /vi. 
Thus, calling 0/v2 the background of the real neutral component of $, the 
tree-level potential is of the form ([5]), with the quartic coupling given by 
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Figure 6: The wall velocity for g = 12 real scalars with \i = 0. 



A = (g 2 + g' 2 ) cos 2 (2/5) /8. Therefore, the tree-level Higgs mass is bounded 
by m 2 H < m 2 z . However, this tree- level relation is spoiled by radiative correc- 
tions (see e.g. [12]) and we shall consider m H as a free parameter. 

In this model, the relevant SM field- dependent masses are those of the 
gauge bosons, 

m 2 w = g 2 <p 2 /A = h 2 w 4> 2 , m| = (g 2 + g' 2 ) <f/A ee h 2 z <p 2 , (55) 
and top quark, 

m 2 (0) = = ^ (56) 

where /it is the Yukawa coupling to H®. We shall work in the limit in which 
the left handed stop is heavy (mq > SOOGeV^). In this case, the one-loop cor- 
rection to the SM is dominated by the right-handed top squark contribution, 
with the field-dependent mass given by 

mf ((f)) ^m 2 u + h 2 ^ 2 , (57) 

where 

/i? = 0.15/4 cos 2/3 + /i 2 [ 1- A ) (58) 

V m Qj 

mfj and rriq are soft breaking parameters, and A t is the stop mixing param- 
eter. If the mass of the right-handed stop is of the order of the top mass 
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or below, the one-loop effective potential f|T2l admits the high-temperature 
expansion 



V (0, T) = D (T 2 - T 2 ) 2 - T IE, 



SMS 



m 



H 



I (8v 2 ) + 5h 2 w /12 + 5/i|/24 + /i?/2 




(59) 



/(4D), 



SM 



where D 

is the cubic-term coefficient in the high-temperature expansion for the SM 
effective potential, E SM ~ (2h 3 w + h 3 z ) /6vr, and M 2 (0) = m| (0) + n t -(T). 
The thermal mass is given by [43J 
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(60) 

where g s is the strong gauge coupling. Following Ref. [15J, we shall set A t = 
for simplicity in the numerical calculation. We shall also take sin 2 /3 = 0.8. 

The phase transition strength is maximized for negative values of the soft 
mass squared mfj ~ —Hf (T) [H], for which the contribution of the term Ai? 
is of the form — E » - - 7 ' •'' 



m lio^l) is ot tne torm —&mssmT4> , with a coefficient E m $sm that may be one 
order of magnitude larger than that of the SM. This would make the phase 
transition sufficiently strong for baryogenesis for Higgs masses as large as 
lOOGeV. However, such large negative values of mfj may induce the presence 
of color breaking minima at zero or finite temperature [15]. Demanding the 
absence of such dangerous minima constrains the Higgs mass to unrealistic 
values. Nevertheless, the two-loop corrections are very important and can 
make the phase transition strongly first-order even for mu ~ [16]. The 
most important two-loop corrections are of the form 2 log and are induced 
by the SM weak gauge bosons, as well as by stop and gluon loops [HH 
In the case of a heavy left-handed stop we have [13] 
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i , Ah 



where the scale depends on the finite corrections and is of order lOOGeV^. 
Following [T3], we will set Ah = lOOGeV^ for the numerical computation, 
given the slight logarithmic dependence of V2 on A#. 

The phase transition takes place at some temperature between the critical 
temperature T c and the temperature at which the barrier between minima 
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disappears. The latter is approximately given by the parameter T in Eq. 
floTI]) . In order to avoid the presence of color-breaking minima, we only con- 
sider values of mfj for which mfj + 11^ (To) > [15J. For the computation 
of the temperature T n we used the high-temperature approximation ( |59i) for 
the one-loop effective potential, together with the two- loop correction ( 16TI) . 
Therefore, for the friction coefficient we also used the high-temperature ap- 
proximation, given by Eqs. © and (jlj). For the stop contribution we have 
rriD ~ mi ~ h{T. 

Figure [7] shows the wall velocity as a function of the stop mass for Higgs 
masses m# = 100, 110, 120, 130 and \4QGeV . In the lower curve (m.# = 
lAOGeV), the cross marks the case n /T„ = 1. To the right of the cross the 
phase transition is weak. As can be seen, for this model the velocity is not 
very sensitive to the stop and Higgs masses. We only obtain deflagrations, 
with velocities in the range v w ~ 0.3 — 0.45. 
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Figure 7: The wall velocity as a function of the stop mass for several Higgs 
masses. The cross in the lower curve marks the point of <p n /T n = 1. 



5.3 TeV fermions 

So far we have considered extensions of the SM for which the relevant con- 
tributions to the effective potential came from new scalars. In Ref. |48j . it 
was shown that in extensions with extra fermions strongly coupled to the 
Higgs field, the phase transition may be sufficiently strong to avoid erasure 
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of the baryon asymmetry in the broken-symmetry phase. Strongly coupled 
fermions, however, make the vacuum unstable unless the Higgs is heavy. In 
the model considered in Ref. [H], this problem was solved by adding heavy 
bosons with similar couplings and number of degrees of freedom, but with 
a large ^-independent mass term, so that they are decoupled from the dy- 
namics at T ~ v. The model can be considered as a particular realization of 
split supersymmetry, where the standard relations between the Yukawa and 
gauge couplings are not fulfilled. Therefore, the fermions are higgsinos and 
gauginos, with a total of 16 d.o.f. 

Depending on the values of the Yukawa couplings and of the mass param- 
eters, the mass eigenvalues can be rather cumbersome. In the simplest case, 
only g = 12 d.o.f. are coupled to the SM Higgs, with degenerate eigenvalues 
of the form m 2 (0) = fi 2 + h 2 (f) 2 . One can assume for simplicity that the 
bosonic stabilizing fields have the same number of d.o.f., and a dispersion 
relation m| (0) = /i| + h 2 <p 2 . We shall also assume, as in Ref. [IE], that 
lis — 0. The maximum value of n$ consistent with stability is obtained by 
requiring the quartic term in Eq. (JH]) to be positive at scales much larger 
than v . Taking into account only the radiative corrections associated with 
the strongly coupled fields, one finds that /i| must be below the value 



In order to minimize the effect of the stabilizing bosons on the strength of 
the phase transition, we will set /is to this maximum value. Notice, however, 
that we have /is ^> /i only for small h, so in general the stabilizing bosons 
are not completely decoupled from the thermal plasma. 

In Fig. El we have plotted the wall velocity for fif = and different values 
of itlh- The parts of the curves on the right of the crosses correspond to 
4> (T n ) /T n > 1. To obtain a strongly first-order phase transition with extra 
fermions, large values of the Yukawa coupling h are needed. Therefore, the 
friction coefficient is in general larger than in models in which the strength 
of the transition is enhanced with extra bosons. As a consequence, the wall 
velocity is smaller. We find velocities v w < 0.25, and as small as v w = 0.05 
for strongly first-order phase transitions. This makes this model interesting 
for baryogenesis, since the generated baryon asymmetry peaks for v w -C 1 
[HE]- Figure [9] shows the case /i ^ 0, for tuh = !25GeV and different values 
of h. In the lower curves, the phase transition is weak on the right of the 
crosses. We see that the wall velocity can either grow or decrease as the 
invariant mass of the extra particles is increased. Nevertheless, v w is still in 
the range 0.05 < v w < 0.25 for strongly first-order transitions. 
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Figure 8: The wall velocity as a function of h, for g = 12 and \x = 0. 
The points to the right of the crosses correspond to phase transitions with 

<f)„/T n > I. 



5.4 Other extensions 

Given the large variety of extensions of the SM, it is not feasible to consider all 
the models. For many models, the strength of the phase transition depends 
only on the particle content. Generally, adding bosons to the SM makes the 
phase transition more strongly first-order, whereas fermions make it weaker. 
To one-loop order, the magnitude of this effect depends on the coupling of the 
extra particles to the Higgs boson, as well as on their number of degrees of 
freedom. The friction also depends on these parameters, increasing with the 
addition of either bosons or fermions. The previous examples correspond to 
such models, in which the extra-SM sector was dominated either by bosons 
or by fermions. For an intermediate model in which bosons and fermions 
have the same values of g, h, and /i, we have obtained intermediate values 
of the velocity, namely, only deflagrations with v w ~ 0.15 — 0.35 for phase 
transitions with <p n /T n > 1. 

On the other hand, the strength of the phase transition may be changed 
without directly changing the particle content of the model. This is what 
happens, e.g., in the MSSM when the two-loop contribution (16T1) is taken 
into account. Another possibility is the introduction of a non-renormalizable 
dimension-six term of the form — v 2 /T\ /A 2 in the Higgs potential, 

which allows to consider a negative quartic coupling |49j . Adding a real 
singlet field S to the SM allows for the possibility of cubic terms of the form 
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Figure 9: The wall velocity as a function of h, for g = 12 and tuh = 125GeV. 
For the points to the right of the crosses the phase transition is weak (0 n /T„ < 

I)- 

(H^H*) S or S 3 in the tree-level potential, which cannot be constructed with 
Higgs doublets. The presence of cubic terms already at zero temperature 
makes it easier to get a strongly first-order electroweak phase transition [36J . 
This possibility exists also in the Next to Minimal Supersymmetric Standard 
Model (NMSSM), which consists of adding a gauge singlet to the MSSM. In 
this model the cubic terms arise as supersymmetry-breaking soft terms. Since 
the strength of the transition is dominated by the cubic terms in the tree- 
level potential, it is not necessary to rely on loop corrections or to consider 
a light stop [50J. 

In these modifications of the SM, the phase transition can be made 
strongly first-order without increasing the friction. In order to explore the 
effect of such a tree-level modification on the wall velocity, we shall consider 
the addition of cubic terms. Considering the full potential makes the model 
considerably more complicated than those we have studied so far, since one 
has to deal with more than one scalar and several free parameters. Instead, 
we shall consider a toy model which consists of adding a term — A<p 3 to the 
tree-level potential (jSJ) for the SM, where A is a free parameter with mass 
dimensions. In this model the parameters of the potential are related to the 
physical Higgs VEV and mass by 2m 2 = Xv 2 — 3Av, m 2 H = 2Xv 2 — 3Av. We 
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shall use the high-temperature potential 

V (0, T) = D SM (T 2 - T 2 ) 4> 2 - (TE SM + A)4> 3 + ^\ (63) 

with the SM values given by D SM = (2h 2 v + 2h\ + h 2 z ) /8, E SM = (2h 3 v + h%) /6tt, 
and Tq = m 2 /D S M- 

We show the results in Fig. [TUJ We have considered values of the pa- 
rameters for which h t <fr/T < 1, so that the high-temperature expansion f[6"3"]) 
is valid. The strength of the phase transition increases very quickly with 
A, and this is reflected in the behavior of the wall velocity. Notice that, 
although in this case the particle content does not change, the friction co- 
efficient increases with the strength of the transition due to its dependence 
on the surface tension. Eventually, the wall velocity stops growing. There is 
also a strong dependence on m#. For smaller values of itlh [pin = 100 and 
200GeV) the strengthening with A is so fast that the deflagration solution 
quickly becomes supersonic and disappears. In these cases, the detonation 
solution does not appear before the phase transition ceases to occur. Inter- 
estingly, detonations appear for weaker phase transitions (higher values of 
m H ). 
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6 Gravitational waves and baryogenesis 



The generation of gravitational waves and of the baryon asymmetry of the 
Universe in the electroweak phase transition seem to be mutually exclusive, 
since baryogenesis requires small wall velocities while GW production re- 
quires large velocities. In this section we discuss which of them is more likely 
in the various extensions of the SM we have considered. 

Both detonations and deflagrations generate gravitational waves [U EJ Qj], 
[T2| \13\ [HI [15j [El [T7J EU [52] . However, higher velocities give stronger signals, 
and Jouguet detonations are often assumed for studying GW generation. As 
shown in Ref. |21j, the Jouguet velocity is not a good approximation and 
detonations, when they exist, are weak detonations. Furthermore, as we have 
seen, in the case of the electroweak phase transition detonations exist only in 
some of the extensions of the SM, and for extreme values of the parameters, 
namely, small Higgs masses and large coupling h. (Interestingly, a non- 
vanishing mass parameter /i may favor the appearance of detonations, as can 
be seen in Fig. H]). In any case, detonations are more likely in models in which 
the strength of the transition is enhanced by tree-level effects, as shown in Fig. 
[TU1 On the other hand, deflagrations with v w > 0.1 may also produce a signal 
of GWs observable by LISA [52] . We remark that the usual assumption that 
detonations are a stronger source of GWs than deflagrations may be wrong. 
Weak detonations have larger velocities than Jouguet detonations, but cause 
a smaller disturbance in the fluid. On the other hand, deflagrations with 
velocities close to c s may cause important perturbations. For most of the 
models we have considered, there is a wide range of parameters for which 
the deflagration velocity is quite large. Extensions of the SM with strongly 
coupled fermions constitute the only exception. 

For the baryon asymmetry of the Universe to be produced at the elec- 
troweak phase transition, the latter must be strong enough, i.e., the condi- 
tion 0/T > 1 must be fulfilled. On the other hand, the generated BAU is 
maximized for velocities v w < 10 -1 , and is strongly suppressed for higher ve- 
locities. This condition demands a weak enough phase transition, how weak 
depending on the friction. These two conditions may restrict significantly 
the parameters of the theory. Models with extra bosons which are strongly 
coupled to the Higgs have been extensively considered, since they easily give 
strongly first-order phase transitions. However, in our examples we have seen 
that these models also tend to give quite large wall velocities. For extensions 
with scalar singlets, we obtained in all the cases v w > 0.2 for 4> n /T n > 1, and 
even larger (v w > 0.5) if tree-level cubic terms are allowed. For the MSSM 
with a light stop we also obtained rather large velocities, v w > 0.35. On 
the other hand, extensions with fermions are not usually taken into account 
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because they tend to weaken the phase transition. However, models with 
fermions can give strong enough phase transitions, as shown in Ref. [IE] for 
the case of fermions with large Yukawa coupling h and heavier stabilizing 
bosons. For this model, with us given by Eq. ( 1621) . we see that the wall 
velocity can be as small ELS Vyj — 5 X KT 2 for <j) n /T n = 1 (see Figs. Eland B>. 

Due to the approximations used in the calculation of the friction coeffi- 
cient, the error in the wall velocity is an 0(1) factor. With a larger friction, 
extensions with bosons may give velocities v w < 0.1 for <p n /T n = 1. In any 
case, it is clear that values of 4> n /T n and v w which are appropriate for baryoge- 
nesis are more easily obtained in models in which the dynamics is dominated 
by fermions. If the friction is smaller than our estimate, we will have higher 
velocities, which favors the generation of GWs. For instance, in the MSSM, if 
we consider a friction coefficient a factor of 3 smaller we obtain deflagrations 
with v w ~ c s and detonations with v w ~ 0.85 for ran = !20GeV. 

7 Conclusions 

We have studied the velocity of bubble walls in the electroweak phase tran- 
sition. We have estimated the friction on the wall due to particles with 
dispersion relation ra 2 = u 2 + h 2 <j) 2 , for arbitrary values of the parameters \i 
and h. We have discussed analytically the cases of small and large /i/T and 
h(j)/T. We have solved numerically the equations for the wall velocity, tak- 
ing into account the friction and the hydrodynamics, and we have computed 
the electroweak wall velocity for several extensions of the Standard Model. 
We have also discussed the implications of our results for baryogenesis and 
gravitational wave generation. 

The friction coefficient, as well as the amount of supercooling, have a 
strong, nontrivial dependence on the coupling h of the particles to the Higgs. 
Furthermore, the friction depends on the bubble wall profile, and thus on the 
strength of the phase transition. Therefore, the wall velocity is not related 
to the strength of the transition in a simple way, and the behavior depends 
on the model. We have found that in general the velocity increases with the 
coupling h, but it can also decrease (see, e.g., Fig. HJ). As we increase the 
mass parameter /i, the friction tends to increase for small /i and to decrease 
for large \i. 

As we have seen, detonations exist only for high values of h and low values 
of run, except in the case in which the strength of the phase transition is due 
to tree-level effects rather than to loop contributions of extra particles. We 
stress that there may be 0(1) errors in our estimate of the friction coefficient. 
For non-relativistic velocities this will not introduce qualitative differences. 
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However, if a model has deflagrations with velocity close to the speed of 
sound, a lower friction may cause the detonation solution to exist and the 
deflagration to disappear. In particular, the interaction rates T need to be 
calculated accurately in each model. 

We have also seen that it is difficult to obtain small velocities in phase 
transitions with <f> n /T n > 1. For instance, the MSSM gives velocities in the 
range 0.35 < v w < 0.45, a hidden sector of scalar singlets gives velocities 
ranging from v w > 0.1 up to supersonic values, and a tree-level cubic term 
causes velocities v w > 0.5. In these models the wall velocity is rather large 
for baryogenesis. On the other hand, such high velocities (either detonations 
or deflagrations) may produce an observable signal of GWs. In contrast, the 
presence of fermions causes smaller velocities, since they increase the friction 
without increasing the phase transition strength. Thus, baryogenesis is more 
likely in models with fermions. In particular, for the model considered in 
Ref. [IE], which contains fermions with large Yukawa couplings h, we have 
found velocities in the range 0.05 < v w < 0.25 for 4> n /T n > 1. 
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